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0(N) linear a model is solvable in the large limit and hence pro- 
vides a useful theoretical laboratory to test various unitarization ap- 
proximations. We find that the large Nc limit and the mo- — > oo limit 
do not commute. In order to get the correct large Nc spectrum one 
^ . has to firstly take the large Nc limit. We argue that the /o(600) meson 

J> \ may not be described as generated dynamically. On the contrary, it is 

^ ' most appropriately described at the same level as the pions, i.e, both 

■ appear explicitly in the effective lagrangian. Actually it is very likely 

! the a meson responsible for the spontaneous chiral symmetry breaking 

\ in a lagrangian with linearly realized chiral symmetry. 

o \ 

^ \ There have been remarkable progresses in recent years in the field of low 

CL|| energy strong interaction physics. In the 1=0, J=0 channel of vrvr scatterings, 

Q^' it is demonstrated that the /o(600) (or a) meson is crucial to adjust chiral 

^ ■ perturbation theory to experiments pj|, and in fact, the /o(600) meson provides 

the dominant contribution to the phase shift at low energies. The Roy type 
equation analyses ^ clearly indicate the existence of a light and broad S matrix 
^ \ pole in the IJ = 00 channel of tttt elastic scattering amplitude, and also in 

' the 1=1/2, J=0 channel nK scattering amplitude {K*{800) or k) [3j. The 

dispersive analyses of Refs. [H O [6] , firmly demonstrate the existence of the a 
and the k resonances as well. AH the results on the pole positions based on 
dispersion techniques agree with each other rather well convincingly, within 
error bars - hence, in the authors' opinion, settling down the long debate 
on whether there exist such light and broad resonances. The a meson can 
also be understood in other approaches. For example, using Dyson-Schwinger 
equations, in Ref. [7] the mass and even the large width of the a meson can 
be determined rather well. 

Having firmly established the existence of sigma and kappa, the next ques- 
tion would be what these resonances are. There exists a long list of references 
with quite different opinions. For example, there are attempts to explain these 
light and broad resonances in linear sigma models at hadron level [8] or at 
quark level [9j. Interestingly, it was observed that poles corresponding to 
these lightest scalar resonances can be found in the unitarized amplitudes cal- 
culated using chiral perturbation theory, i.e., not having to include the scalars 
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explicitly in the effective lagrangian. Hence the lightest scalars are considered 
as 'dynamically generated' [10]. Unlike a typical resonance, for example, the 
p(770) meson, it is found that the a and k resonance poles move to infinity 
on the complex s plane in the large A^^^ limit [11]. This different N^. behavior 
is hence believed to be an evidence to support that the a and k are generated 
dynamically. This aim of this paper is to examine the latter approach closely 
and try to clarify the meaning of 'dynamically generated' resonance. 

The idea that a particle can be generated dynamically is not new. Basde- 
vant and Lee have succeeded in generating a dynamical p resonance from the 
renormalizable linear sigma model [12] . On the other side, Basdevant and Zinn- 
Justin studied the massive Yang-Mills lagrangian in which the p meson acts 
as the massive gauge boson. Through a proper unitarization procedure, they 
generate in d wave the /2(1270) resonance and a broad s wave a resonance [T3] . 
These activities are closely connected to the old concept of bootstrap duality: 
in one lagrangian model, a particle may be called elementary (i.e., exhibits 
itself explicitly in the lagrangian), others are called dynamical (i.e., only show 
up as 5* matrix poles but not included in the lagrangian); whereas in other 
lagrangian models, the particle's role as 'elementary' or 'dynamical' can be 
just reversed. 

The idea of bootstrap duality is very appealing. However it is difficult 
to justify it in practice, since it involves nonperturbative calculations of field 
theory. The only way at present is to use unitarization approximations to 
investigate it. One of the most frequently used method is the Pade approx- 
imation. Nevertheless it is difficult to justify to what extent one can trust 
the predictions of Pade approximation, especially when it is applied to chiral 
perturbation theory (xPT) amplitudes [Mj. Concerning the difficulty of the 
problem, in this letter we will firstly focus on the 0{N) linear sigma model. 
This model is solvable in the large N limit [15] and hence provides a good the- 
oretical laboratory to test the validity of various unitarization approximations. 
From these studies, one can gain better understandings to the unitarization 
approximation, when it may be trusted and to what extent they can be trusted. 

The lagrangian of linear 0(N) model is. 



where $ = ($i, $2, ■ ■ ■ , ^n)'^- We work in the symmetry breaking phase, so 
< and the $ field develops a VEV, v"^ =< >= -2Niv? /\ which 

is proportional to N (i.e., number of flavors). To be more realistic, we discuss 
the TTTT scattering in 0(N) model with massive pions, by adding an explicit 
symmetry breaking lagrangian to the above £, 
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jO^S.B. = C^N 



(2) 



and c = vm. 



4. The pseudo-goldstone boson scattering amplitude is. 



(3) 
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where 



and also 

^ob', ml) = -2 / (2vr)i5 (p _ ^Ddk + pf - ml) ' 
The integral in above is made finite using dimensional regularization, 

f d^k 11 i ~ 

I (27r)^P-m2 [p + k)^-ml ^ ~ 16^^°^^' ""'^ ' 

with {e = 2-D/2) 

A = -i^[r(e)+ln47r + ln4] , (7) 
lovr^ mi 



I ^uit — s „ 1 



4ml- 



B,{s, mi) = -2 - «\/ ^7-753= + 



In order to compare with the calculations made later in this paper, notice 
that here we adopt a different strategy to the 0{N) model: we treat it as a 
cutoff effective theory, which means all the parameters in lagrangian are taken 
as finite. In cutoff effective theory one introduces a cutoff A in the above 
divergent integral. This can also be done using dimensional regularization 
scheme, by introducing the cutoff A in the following manner, 

r(e)+ln47r + ln^ ^In— . (9) 

mi mz 

TV 7T 

One can demonstrate that the above procedure is physically equivalent to 
using a sharp momentum cutoff method, when the energy scale is well below 
A. Notice that, as clearly shown from Eqs. (jlj) and ([5]) that the notorious 
square divergence does not occur in the scattering amplitude. The reason is 
that the square divergence only occurs in the effective potential, and once it 
is absorbed by the bare a mass through mass renormalization it disappears 
entirely everywhere else. This observation justifies the negligence of square 
divergence in the latter calculation in effective theory. Another issue with 
respect to the 0{N) model is that the 0{N) model may contain a tachyon in 
some regularization schemes. The tachyon also appears in the present cutoff 
version. But it is numerically verified that the tachyon pole locates at a scale 
much larger than the cutoff itself and is hence of no concern. In the cutoff 
version the parameter A appeared in Eq. (jlj) is the same as the bare parameter 
A appeared in the lagrangian Eq. ([1]) and obeys the relation 

,_ ml- ml 
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The denominator of the propagator D is given by the following expression, 



1 1 771^ 1]^ 

G^is- ml){- - ^,{Bois,ml) + In ^]} - - . (11) 

Poles appeared in vrvr scattering are determined by solving the equation G^^ = 

0, where G^^ is the analytic continuation of G on the second sheet. From above 
we realize that the true sigma pole position is a function of the bare sigma mass 
parameter, for fixed cutoff A. The sigma pole trajectory is plotted in Fig. 1. 
Throughout this paper we take v/\/N = = 93MeV and A = 1.5GeV. The 
bare mass of sigma in Fig. 1 ranges from 400MeV to 1200MeV, and each value 
differs by lOOMeV. 

In the following we will discuss variations of the above 0(N) lagrangian 
and calculate the scattering amplitude using different unitarization approxi- 
mations. By comparing with the above rigorous results, we can gain some 
knowledge on the quality of various approximation methods. It is worth notic- 
ing that, for the toy linear 0(N) model, it is demonstrated that the [n,n] Fade 
amplitudes gives the exact sigma pole location, as comparing with the exact 
solution [16j . The K matrix approach does not reproduce the exact sigma pole 
location, but it is still a good approximation [16]. In the following we will show 
that the nice property of these approximation methods will not be maintained 
if the pion fields are expressed in the non-linear representation. 

The next thing we will do is to make polar decomposition of the 0(N) field, 

1. e., 

$ = y?0, 0= (01,02, 07v) (12) 

with constraint 

0? + 02 + 03 + + 0^ = 1 • (13) 

In the symmetry breaking phase we define the shifted field, 

!f = a + v (14) 

and normalize the field as 

vf = t>0 . (15) 
With above preparations we can obtain the polar decomposition lagrangian. 



Cpoiar = C^ + ^il + -)'(9^vr ■ d'^fr + d^Vv^ - 7f ■ WVv^ - vf ■ 7f) (16) 



with 



In the above formula we have used the constraint condition, Eq. f[T3|) . so vf 
{tti,tt2, ttn-i)- Now the symmetry breaking lagrangian, Eq. ([2]) reads. 



= Id.ad^a - \m\a + vf - ^{a + vf . (17) 



i^SB = c{a + v)\ 1- 



TT ■ TT 

y2 
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Then we expand the square root to get the interacting vertices which are 
needed for evaluating O(p^) vrvr tctc processes. 



2 V t>^ — TT-TT Vf 

2 V 2v^ 2v^ 

H jfTT ■ (9uVf)^7r • TT ^TT • TT ^(TTT ■ TT ^{tI ■ TxY ^(vf ■ vf)'^ 

2v^^ ^ ' 2 2v Sv^^ ' 16v^^ ' 

+ ■■■. (19) 

Notice that in the above lagrangian the a^vr^ term does not contribute to tttt 
scattering amphtude at 1-loop level since its contribution is 1/N suppressed. 
There are other variants of the above effective lagrangian. For example we 
may simply neglect the a field in Eq. f|T6|) and get the 0(N) non-linear sigma 
model 

1 



i^O{N) non-linear = 2^^^^ ' "^^^ + <9/. V v"^ - TT ■ 7x0'' V - TT • vf) . (20) 

Expanding the pion fields from the above lagrangian one gets (also obtainable 
from Eq. (IT^ by simply neglecting the sigma field): 



-NL = -f^M^ 



2 2 2i;2^ ^ ^ 8v^^ ' 



+ ^(vr.9,vr)^vr.vr-^(7r.vr)^ + ... . (21) 

Also one can integrate out the a field at tree level from Eq. fll6p to get a 
modified 0(N) non-linear sigma model: 

This lagrangian resembles the O(p^) chiral perturbation theory lagrangian. 

Starting from Eq. ( fT9i) we can calculate the HaTXt — * vTcVTrf scattering ampli- 
tude up to O(p^) in large N limit: 

1. The contribution from tree level 47r contact term, or current algebra 
term, denoted as T2 (~0(p^)), 

T2 = SabScdA2{s) + SacSbdA2(t) + (5ad5{,c^2 («) (23) 

where 



A2{s) = -^{s - ml) . (24) 
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2. The contribution from tree level O(p^) term contained in Eq. (122|) . de- 
noted as T4, is 

74 = SabScdA^is) + SacSbdA^it) + SadSbcA4,{u) (25) 

where 

^4(.) = ^-^^^ . (26) 

3. The contribution from tree level a exchange in Eq. f|T9|) (~ 0{p'^)), 

Ta = 6abScdAa{s) + SaJbdAait) + 6adSbcAa{u) (27) 

where 

Ms) = 4(1^ . (28) 

which is O(p^) at low energies. 

4. The contribution from 1-loop diagrams (in large N limit and including 
only O(p^) terms): 

Tloop = Sab^cdAioopis) + SacSbdAioopit) + Sad^bcA-loopiu) (29) 

where 

A^.. 0.0 „ , 3m! ~2s r d^k 1 



Aioop = ^[(s - ^l) Bo{s- ml) - j 



(27r)^A;2 -m2J ' 
(30) 



The quadratic divergent integral from bubble diagram in above is treated in 
the standard way in dimensional regularization scheme0 

"^""^ ^ im^A - -^Aim^) , (31) 



(27r)^ P - m2 167r2 
where 

A{m^) = -m^ . (32) 

As before we introduce the cutoff in the same way as in calculating the solvable 
linear 0{N) model, i.e., 

r(e) + ln47r + In — ^ In — . 



^If we were to use a square divergence, i.e., a to replace the following integral, the 
resulting pole location will be disastrous as comparing with the original 0{N) model pre- 
diction. Also a naive cutoff regularization scheme spoils chiral symmetry and hence is not 
correct. An effort to remedy the deficiency of the naive cutoff regularization is discussed in 
Ref. [m. 
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One can check that the above treatment to ignore the square cutoff leads to 
the same result in the chiral limit to 1-loop amplitude comparing with the one 
obtained from the original 0{N) model. 
Now consider the 1=0 channel, 

T^=°(s, t, u) = {N- l)A{s) + A{t) + A{u) . (33) 

In the large N limit, only the s channel amplitude A{s) term survives. The 
J=0 partial wave is simply 

TOO ^ 1 iVA(s) . (34) 

OZTT 

With the above preparation we will in the following construct various uni- 
tarized amplitudes based on tree and one loop calculations using different 
lagrangians. The unitary amplitudes being constructed are: Ki (using tree 
level amplitude in pure 0{N) non-linear a model), K2 {Ki plus a exchange), 
K3 (1-loop calculation in the pure 0{N) non-linear a model), K4 (1-loop with 
a exchange included using polar decomposed lagrangian); as well as Pi ([1,1] 
Pade for the pure 0{N) non-linear a model), P2 ([1,1] Pade with the a being 
integrated out, i.e., of Eq. (122|) ). P3 ([1,1] Pade with explicit a exchange, i.e., of 
Eq. (122|) ). For simplifying the expressions we in the following use v (=93MeV) 
to replace v/^/N. In the following we will always work in the IJ=00 channel 
and hence for simplicity we drop out the superscripts 00 in the partial wave 
amplitude. 

Ki amplitude: This approximation corresponds to the simplest K matrix 
amplitude in the pure 0(N) non-linear sigma model. Also it is the simplest K 
matrix amplitude one obtains from Eq. f|T9l) . One has 

1+tpK 



Here Ki = T2 and hence 



1 + ipT2 , . 

S = , 36 

1 - ipT2 



where p = Jl — Ami I ^ ^'^d 



T2 = ,^^{s-ml). (37) 

The djTiamically generated second sheet poles correspond to the zeros of the 
5* matrix, we show the pole position in Fig. 1 labeled as Ki. Notice that 
there is no free parameter here. We see from Fig. 1 that Ki is not a good 
approximation numerically to the real solution (i.e., the solution of the linear 
sigma model). There is another serious disadvantage in this approximation, 
which can be clearly seen in the chiral limit. That is the pole location, is 
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proportional to v and hence is proportional to the square root of the number 
of colors, Nr,: 



^ oc V^327r/2 ^ ^jv^ _ (38) 

So when Nc ^ oo this pole disappears! On the other side the sigma pole in 
the linear sigma model behaves as mo- ~ 0{N^) in power counting. 

K2 amplitude: K2 = T2 + T^. This approximation corresponds to calculat- 
ing the scattering amplitude at tree level in the polar decomposition lagrangian 
Eq. (fT9|) . and it is also equivalent to the lowest order K matrix amplitude eval- 
uated in the linear sigma model|§ We have, 

with 

= ^^^^ . (40) 

d/TTf s — mi 



and 

K2 = T2 + T, 



11 1 
ml ml — ml 



327rf 2 s — m^ 

The second sheet poles correspond to the solution of the following equation. 



/ s — 4?n^ X 

s-ml + iJ ^-^(^rnl-ml){s-ml) = . (41) 

V s oZttv 

The pole trajectory is plotted in Fig. 1 labeled as K2. Apparently K2 works 
much better than Ki. Especially the pole position is now ~ 0{N^). But 
numerically it is still a rather poor approximation as can be seen from Fig. 1. 

By analyzing the K2 approximant one can learn an important lesson. It is 
easy to find solution of Eq. fHTl) in the chiral limit: 

"^P°'e = j.n,2 • (42) 

327^)2 

In the above equation if one formally sends m„ — >■ 00 one recovers the solution 
Eq. (l38ll in the chiral limit, i.e., Spoie = i32nv'^ — ^ 00 when — >■ 00. Neverthe- 
less if one firstly take the large limit in Eq. P2l) one gets simply Spoie = rnl. 
Apparently the latter one is the correct relation we are looking for. Therefore 
we conclude that one must take the large Nc limit before integrating out the 
sigma degree of freedom in order to get the correct large spectrum. 



^In the linear sigma model one expands the amplitude according to the power of the 
coupling constant A and constructs a K matrix amplitude using the tree level perturbative 
amplitude. 
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amplitude: = T2 + Re[T;oop]. This approximation corresponds to a 
1-loop calculation in the pure 0(N) non-hnear sigma model. Notice that here 
we adopt the cutoff effective lagrangian approach, therefore no renormalization 
procedure and hence no counter terms are needed. 

l-ip{T2 + MTloop]) ^ ' 

with 

+(3m^-2sK(l + ln^)}. (44) 

Notice that in this approximation the sigma pole is dynamically generated and 
its location is also fixed. The pole position is shown in Fig. 1 labeled as K^. 
Notice that here, hke the Ki case, one also has ~ 0{N^/'^). 

amplitude: = T2 + T^. + Re{Tioop)- This approximation corresponds 
to considering all contributions up to 0{p^) term in the polar decomposed 
0(N) model. We have 

1 - ip[T2 +T^ + Re{Tioop)] ■ ^ ' 

The pole trajectory is plotted in Fig. 1, labeled as K4. Comparing with K2, 
the quality of the result is not actually improved. What is even worse here is 
that the amplitude contains a spurious pole on the second sheet, not far 
from the sigma pole position. 

We conclude in general that the K matrix unitarization approach to the 
amplitude obtained from a lagrangian with non-linearly realized chiral sym- 
metry docs not work very well. Still, there exist other distant spurious poles 
(can be on both sheets) in these K matrix amplitudes, which are not shown in 
Fig. 1. On the contrary, the K matrix approach works fairly well on the linear 
0(N) model. 

In the following we turn to discuss the Fade approximations which are very 
frequently used in phenomenological discussions. Notice that Fade amplitudes 
also preserve one bad heritage of the original exact amplitude - the tachyon 
pole. But as stated before the appearance of the tachyon pole does not matter 
in practice in the phenomenological discussions. 

Pi amplitude: This approximation corresponds to the [1,1] Fade approxi- 
mation in the cutoff pure non- linear 0(N) model, 

T- rr. % ■ (46) 

— J- loop 
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The pole position is shown in Fig. 2 labeled as Pi. Certainly it is better than 
Ki but similar in quality with K3. The same problem remains, however, that 

is ^ (XV (X 0{Nl/^). 

P2 amplitude: This approximant corresponds to a 1-loop calculation using 
/:^of Eq. m- 

T = ? (47) 

2 — J-i — loop 



where 



(s - ml^'^ 



327rf 

Second sheet poles are determined by solving the following equation 



s-mi- ^ - 77:-^{(s - mlf[2 + In — + p(ln — ^ - zvr) 

A2 



+ {3ml-2s)ml{l + ln^)} = (49) 

The pole trajectory is shown in Fig. 2 labeled as P2- This approximation is 
the best among all. It is not difficult to examine using Eq. ( l49l) that if one 
takes the large Nc limit one finds ~ rria. On the other sides, if one sends 
rrifj 00 before taking the large Nc limit one gets oc f ~ O^N^^'^). Once 
again we find that the mo- 00 limit and the large Nc limit do not commute. 



P3 amplitude: This approximation corresponds to the [1,1] Fade approxi- 
mation up to O(p^) term in the polar decomposition lagrangian Eq. (fT9|) . We 
have, 

rp2 

1 rji rji ' (^^} 

2 — J-a — J- loop 

The pole trajectory is plotted in Fig. 2 labeled as P3. Notice that adding a 
sigma pole explicitly in the lagrangian does not improve the quality of output. 
What is even worse is that it also predicts a spurious pole not very far from 
the sigma pole. Our analysis seems to suggest that a polar decomposed sigma 
model lagrangian works even worse than the situation when the sigma field 
is completely integrated out. It is however not clear whether this observation 
remain valid in the realistic case. 

The pole appeared in the unitarized amplitude is dynamical in the case of 
0{N) non- linear a model, and also dynamical in the case of Eq. fl22|) . In the 
former case the pole moves to infinity as A^^c — + 00 whereas in the latter case 
the pole falls down to the real axis when Nc —>■ 00. For the latter we know 
that the dynamical pole is nothing but just represents the original a pole being 
integrated out from Eq. (fT9!) or Eq. ([1]), as is revealed by the simple exercise 
on P2 amplitude, hence 'elementary' in its origin. The drastically different Nc 
dependence of the pole trajectories, in our understanding, merely reflects the 
different order when taking the large Nc and m„ 00 limits. 
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In a more realistic calculation however, the inverse problem becomes prob- 
lematic since the Pade approximants for chiral perturbation theory bring, to- 
gether with the sigma pole, also spurious poles and some of which violate 
analyticity [H]. Nevertheless the toy model we discussed here well resembles 
many aspects of the real situation. For example we pointed out that in the 
simplest Ki approximant, the dynamical pole being generated has the deferent 
Nc behavior comparing with the a meson in the linear 0(N) model. The situ- 
ation is very similar in realistic case. In Ref. [1^, based upon a [0,1] Pade or 
Ki approximant the authors suggest that the current algebra term dominates 
at low energy, since the pole position of the sigma meson obtained as such is. 



and is found to be not far from the pole position obtained from a more careful 
numerical analysis. According to the analysis made in this paper, even though 
the current algebra term dominates low energy physics, the Nc dependence of 
the pole trajectory in Eq. flHT]) does not in any sense prove that the a pole is 
dynamical. The current algebra term in Eq. (fTOjl l^ gives exactly the same Nc 
dependence as in Eq. f lSlll i meanwhile it comes from the linear sigma model 
with explicit a degree of of freedom! If Eq. (15T|) were correct for arbitrarily 
large Nc, the pole will move to infinity when Nc oo, hence contradicting 
the Nc counting rule for a conventional meson state. If this were true the 
sigma meson had to be a dynamical pole. Going beyond the leading current 
algebra amplitude, contrary to what is claimed in Ref. [19], in the [1,1] Pade 
approximant in reality, the a pole actually finally falls on the real axis. That 
is the pole position as given in Eq. (IHTj) receive important corrections when Nc 
is large: 



valid in the large Nc and chiral limit. It is worth noticing that it was not 
clear what is the approximation made to obtain Eq. ( l52l) . Using the PKU 
parametrization form [H El E] one explicitly demonstrates that the Eq. ( l52i) is 
obtained under two further assumptions: 1) assuming one pole (the a pole) 
dominates in the s channel at low energies; 2) neglecting all the crossed channel 
resonance exchanges [20]. It is understood that the assumption 1 actually 
fixes the Nc behavior of the a pole, i.e., rria- ~ 0(1), To- ~ 0{1/Nc) [SU 
[20] . The assumption 2 does no harm in distorting the large Nc behavior. 
However, these assumptions are only approximately correct when A^^^ = 3 and 
fail at large Nc, as demonstrated in Ref. [22]. Actually it is likely that at 
large Nc, crossed channel resonance exchange dominates low energy physics. 

■^Or more precisely, when apply Ki approximation to Eq. (|19p . 

"^There is a factor of 2 difference between the coefficient in Eq. ([STjl and Eq. ([38]) . 




(51) 



1 + 16m/2A ' 
A = :^{22L, + UL2 + IIL,) oc 0{N'c) , 



(52) 
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Hence Eq. (152|) can neither be used to discuss the large Nc behavior of the 
a pole trajectory. In Ref. [19], it is observed that, unlike the p meson which 
almost falls straightforwardly down to the real axis when Nc increases, the a 
meson is very unwilling to fall down. Based on this observation the authors of 
Ref. [IS] concluded that the a pole they found from their unitarized amplitude 
a dynamically generated one. Even if Eq. fl52p is invalid to be useful when 
discussing the large behavior of the a, the discussions given in this paper 
make it clear that the bent structure of the a pole trajectory should not be 
used to argue that the a is dynamical. We plot in Fig. [3] the Nc dependence 
of the a pole in linear 0{N) model (as well as the P2 approximant). Fig. [3] 
clearly shows that, even the 'elementary' a does not fall down to the real axis 
straightforwardly. Though not a rigorous proof, in our opinion the similar 
bent structure of the a pole trajectory both in Fig. [3] and in the unitarized 
chiral perturbation theory amplitude suggest one simple thing: the observed 
/o(600) state is the a meson being responsible for the spontaneous breaking 
of chiral symmetry. In 1+1 dimensions, using an exact solution for the large- 
N limit, Ref. [23] showed dynamical generation of the sigma meson starting 
with the nonlinear sigma model. In 4 dimensions the situation remains to 
be understood. In any case in the linear a model the a meson is the chiral 
partner of the pseudo-goldstone bosons and the two mix each other under chiral 
rotations. It suggests that the a meson shares one dynamical property of the 
pion, i.e., being a collective mode. Hence the appropriate way to describe it is 
not to consider it as dynamically generated from pion fields, rather it should 
be discussed at the equal level as pions in the effective lagrangian approach. 
Also the a meson should not be properly described by any simple quark models 
since a pion can not be. 

This work is supported in part by National Natural Science Foundation of 
China under contract number 10575002 and 10421503. 
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Figure 1: The a pole position from different K matrix unitarizations comparing 
with the rigorous solution. The open square denotes the exact position of the 
0{N) model, others are explained in the text. 
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Figure 2: The a pole position from different Fade approximations comparing 
with the rigorous solution. The open square denotes the exact position of the 
0{N) model, others are explained in the text. 
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Figure 3: dependence of the a pole position. Full circle: linear 0{N) sigma 
model; open circle: P2 approximation. In drawing the plot N^. starts from 3 
and each pole is differed by increasing A^c by 1. The bare mass m„ is taken as 
1.2GeV. 
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